Scaling behavior of fragment shapes 
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We present an experimental and theoretical study of the shape of fragments generated by explosive 
and impact loading of closed shells. Based on high speed imaging, we have determined the frag- 
mentation mechanism of shells. Experiments have shown that the fragments vary from completely 
isotropic to highly anisotropic elongated shapes, depending on the microscopic cracking mechanism 
of the shell-material. Anisotropic fragments proved to have self-affine character described by a scal- 
ing exponent. The distribution of fragment shapes exhibits a power law decay. The robustness 
of the scaling laws is illustrated by a stochastic hierarchical model of fragmentation. Our results 
provide a possible improvement of the representation of fragment shapes in models of space debris. 

PACS numbers: 46.50.+a, 62.20.Mk, 64.60.-i 



Spacecrafts and satellites, during their mission and ser- 
vice time, are exposed to the danger of impact with pieces 
of space debris, which is a growing population of rocket 
bodies, non-functioning spacecrafts, rocket fuel ejecta 
and pieces of fragmented material accumulated during 
40 years of space exploration pj. In order to minimize 
the potential hazard, objects of size larger than 10 cm are 
continuously tracked in space and their orbits are taken 
into account for space activities. Fragmentation events 
like on-orbit explosions of fuel containers of upper rocket 
stages and secondary breakups of fragments due to mu- 
tual collisions, are the main source of the proliferation 
of space debris, creating a large number of small frag- 
ments which cannot be tracked. For safety reasons it is 
essential to work out models of fragmentation, i.e. the 
breaking of objects into smaller pieces, which are able to 
predict the consequences of on-orbit explosions and im- 
pacts @. The NASA breakup model EVOLVE05 0, also 
implemented by other space agencies, represents the frag- 
ments in terms of their characteristic length L c , surface- 
to-mass ratio A/m and velocity v. Other quantities like 
the fragment mass m are determined from scaling rela- 
tions. Model calculations are performed in a phenomeno- 
logical way, i.e. based on experiments and on-orbit ob- 
servations of breakup events, probability distributions of 
the above quantities are prescribed. Monte Carlo sim- 
ulations are carried out taking into account the specific 
initial conditions of the event studied 0] ■ The orbits of 
fragments are determined by their velocity, however, the 
lifetime of the orbits is limited by the atmospheric drag 
which mainly depends on the shape of the fragments. 
The probability of impact of debris pieces with a space- 
craft and the resulting damage can be calculated from 
their velocity, mass, size, and shape. The precision and 
predictive power of model calculations strongly rely on 
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FIG. 1: Consecutive snapshots of an exploding hen egg-shell 
at a frame-rate 15000 1/sec. Merging of side- branches of 
the hierarchical crack-tree and cracks leading to secondary 
breakup are indicated by arrows in b) and d), respectively. 

the quality of the input distributions and the validity of 
scaling relations used. 

General studies on fragmentation phenomena mostly 
focused on the understanding of the mass distribution 
F(m) of fragments. For bulk solids, power law fragment 
mass distributions F(m) ~ mT T have been obtained un- 
der widely varying conditions with universal exponents r 
depending mainly on the spatial dimension d BL IH • 
Recently, we have pointed out that the fragmentation of 
shell-like objects, like fuel containers or rocket bodies rel- 
evant for space debris, forms an independent universality 
class 0, , which is also supported by other studies 0] . 

In this paper we present a study of the shape of 
fragments generated by explosive and impact loading of 
closed shells. High speed imaging investigation of the ex- 
plosion process showed that the shape of fragments is de- 
termined by the underlying cracking mechanism of shells 
which strongly depends on microscopic material proper- 
ties: fragments of isotropic shape are obtained for mate- 
rials where the branching-merging mechanism of cracks 
governs the breakup, while the formation of long straight 
cracks results in fragments of a high degree of anisotropy. 
Fragments of anisotropic shape proved to have sclf-affinc 
character described by a scaling exponent. The distribu- 
tion of fragment shapes has a power law decay with a ma- 
terial dependent exponent. We illustrate the robustness 
of the scaling laws of the shape of shell fragments by a 
hierarchical stochastic model. Our results suggest a pos- 
sible improvement of the description of fragment shapes 
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FIG. 2: Consecutive snapshots of an exploding hollow glass 
sphere at a frame-rate 15000 1/sec. Hot spots can be seen in 
a) radiating a large number of straight cracks b, c, d. 

in phenomenological breakup models of space debris by 
clearly separating bulk and shell fragmentation and by 
taking into account the effect of the cracking mechanism 
of different shell materials on fragment shapes. 

In order to understand how the breakup mechanism of 
shells determines the shape of fragments, we have carried 
out explosion experiments of closed shells made of brittle 
materials with different microscopic properties such as 
hen and quail egg-shells and hollow glass spheres. The 
egg-shells were used as a cheap and easy-to-handle brit- 
tle bio-ceramics with a highly disordered microstructure 
[Toj . A Photron APX ultima high speed camera with 
frame rate 15000/s and spatial resolution 256 x 256 pix- 
els was used to follow the time evolution of the explosion 
process, which also enabled us to study the dynamics of 
crack formation on the surface, and for the first time in 
the literature provided direct access to the mechanism 
of fragmentation. The shells were filled with a stoichio- 
metric hydrogen-oxygen mixture which was electrically 
ignited approximately in the center of the shell. The 
analysis of the explosion of 20 egg-shells showed that the 
breakup process starts with the nucleation of a few cracks 
at the flatter end of the egg (see Fig.^l. Since the energy 
stored in the expanded shell at the instant of crack nucle- 
ation is high compared to the energy released by the free 
fracture surface, the cracks propagate at a rather high 
speed. The instability of the propagating crack results 
in sequential splitting of the crack tip at almost regular 
distances, trig gered by the heterogeneities of the material 
(Fig. fill ]. The propagating sub-branches accelerate 
due to the expansion of the shell and can undergo fur- 
ther splittings giving rise to a hierarchical tree-like crack 
pattern (Fig. \\%>) . Fragments arc formed along the main 
cracks by the merging of adjacent side branches at al- 
most right angles (Fig. ^). At most 3 — 4 splittings 
can be observed along a main branch, however, merg- 
ing typically occurs at the first two levels of the hierar- 
chy. We note that the area of fragments generated at 
this stage does not show large variations, it is practi- 
cally determined by the inherent length scale of crack tip 
splitting, furthermore, the shape of fragments is more 
or less isotropic. The branching-merging process is ini- 
tially governed by the in-plane deformation of the shell, 
however, as the expansion increases the out-of-plane de- 
formation dominates giving rise to further cracks mainly 
perpendicular to the former ones (Fig.^i). This cracking 
proceeds again in a sequential manner, typically breaking 
the fragments into two pieces until a stable configuration 



is reached. 

The cracking mechanism discussed above should be 
generic to materials with a strongly disordered mi- 
crostructure Q. However, for shells made of materials 
like glass.the breaking mechanism can be significantly 
different 0, ■ Experiments on exploding hollow glass 
spheres have shown that the breakup starts at one hot 
spot with random position on the surface (Fig. ^i) from 
which long straight cracks radiate without any appar- 
ent branching (Fig. b). This cracking mechanism re- 
sults in a large number of long thin fragments having 
also a relatively large curvature which makes them un- 
stable against bending. During the expansion of the 
sphere these primary fragments undergo a sequential 
breakup process due to the out-of-plane bending defor- 
mation (Fig. [St, d). Further details on the dynamics of 
cracking of shells will be provided in Ref . 01 ■ Egg-shells 
and hollow glass spheres were also fragmented by impact 
with a hard wall, which produced the same type of frag- 
ments @- 

In the final state of the breakup process, the frag- 
ments were carefully collected and digitized with a scan- 
ner for further evaluation. It can be observed in Figs. ^ 
121 and in the inset of Fig. that for the different types 
of materials considered, the fragments are always com- 
pact two-dimensional objects with little surface rough- 
ness, however, their overall shape can vary from com- 
pletely isotropic (egg-shell) to highly anisotropic (glass) 
depending on the cracking mechanisms. The mass m 
and surface A of fragments is defined as the number 
of pixels N and the contour length of the spots in the 
digital image, respectively. We characterize the linear 
extension of fragments by their radius of gyration as 
Rg = (l/N)J2i^j=i {^i — rj) 2 , where the sum goes over 
the TV pixels fi of the fragments. In order to reveal how 
the shape of fragments varies with their size, in Fig. |3| 
the average fragment mass (m) is presented as a function 
of R g for different materials from impact and explosion 
experiments. It is important to note that in all cases 
power law functional forms (771) ~ Rg are obtained with 
a high quality, however, the exponent a depends on the 
structure of the crack pattern. Since the egg-shell pieces 
have regular isotropic shape, their mass increases with 
the square of R g and hence a = 2 ± 0.05 was fitted. The 
large glass fragments are characterized by a significantly 
lower value of the exponent a — 1.5 ± 0.08, while for 
small glass pieces one observes a crossover to isotropic 
shape with a = 2 ± 0.08. The value a < 2 implies that 
the fragments have self-affine character, i.e. the larger 
they are, the more elongated they get. Note that similar 
anisotropy and self-affinity of fragment shapes was not 
observed in d dimensional bulk fragmentation (d = 2, 3) 

sin. 

To quantify this behavior, let us consider that the shell 
fragments have a rectangular shape with side lengths 
a and b, hence, the surface, mass, and radius of gyra- 
tion can be obtained as A = 2(a + b), m = ab, and 
R g = \J a 2 + 6 2 /(2-\/3), respectively. The fragment mass 
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FIG. 3: Log- log plot of (m) as a function of R g . The arrow in- 
dicates the crossover point of glass fragments. Inset: scanned 
pieces of egg (left) and glass (right). 
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FIG. 4: Distribution / of the shape parameter S. Inset: dis- 
tribution g of A/m. The legend is the same as in Fig. El 



can be expressed in terms of the aspect ratio r = a/b 
and R g as m ~ R 2 /(r + l/r), where even for moder- 
ately elongated fragments the approximation m ^ R 2 /r 
is valid. Consequently, for fragments with exponents 
a < 2, the aspect ratio must increase as a power of R g so 
that r ~ R s g . Hence, m ~ Rg~ S follows, and a = 2 — S, 
where 5 w 1/2 was obtained in the experiments. It is in- 
teresting to note that the value 5=1/2 has been found in 
a broad class of systems producing self-affine structures, 
for instance, for the scaling of the width with the length 
of the arms of noise reduced DLA clusters, for clusters of 
directed percolation or for the Holder exponent of one- 
dimensional random walks |15| . 

Apart from the shape of individual fragments, it is also 
important to know the probability of occurrence of a spe- 
cific fragment shape in the final state of a breakup pro- 
cess. The NASA breakup model characterizes the shape 
of fragments by the surface-to-mass ratio A/m, the dis- 
tribution of which is fitted by a linear combination of 
Gaussian distributions The functional form of the 
corresponding distributions g(A/m) of our shell pieces in 
the inset of Fig. 0] again shows a strong dependence on 
the cracking mechanism. For isotropic fragments a rea- 
sonable fit could be obtained with Gaussians in agree- 
ment with the NASA model Q, however, for anisotropic 
fragments g increases monotonically. The small sized 
isotropic and the large very elongated anisotropic frag- 
ments both have large A/m value which prevents clear 
shape identification. To obtain a better characteriza- 
tion of fragment shapes, we introduce a dimensionless 
shape parameter 5 defined as S = ^R g , multiplying 
the surface-to-mass ratio A/m by the radius of gyration 
R g . Assuming rectangular objects, the shape parameter 
S takes the form S = (a + b)Va 2 + b 2 / (V3ab) . For frag- 
ments of isotropic shape a « b, it follows that S ~ 1.63, 
which is indicated by the vertical dashed line in Fig. 0] If 
the fragments are elongated a>i, the shape parameter 



S ~ a/b coincides with the aspect ratio r characterizing 
the degree of anisotropy. Corresponding to the crack- 
ing mechanisms, the distributions f(S) of different ma- 
terials and fragmentation modes (explosion and impact) 
form two groups in Fig. Fragments of a low degree 
of anisotropy, irrespective of their size, contribute to the 
maximum of f(S) in the vicinity of S ~ 1.63. Since 
egg-shell fragments are mostly isotropic at all sizes, the 
distribution f(S) decreases rapidly over a narrow inter- 
val of S. The remarkable feature of the distribution / 
is that it follows a power law decay f(S) ~ S^ 13 , where 
the exponent /3 = 6.8 ± 0.3 is obtained for isotropic frag- 
ments over a limited scaling range, while (3 = 3.5 ± 0.2 
follows when the cracking mechanism favors the forma- 
tion of anisotropic fragments. 

Discrete element models of shell fragmentation usu- 
ally consider highly disordered brittle materials and pro- 
vide isotropic fragment shapes, but they have difficul- 
ties to capture microscopic mechanisms resulting in long 
straight cracks 0, El- We propose a simple stochastic 
binary breakup model in the spirit of Refs. 0, 0] to 
better understand the experimental findings. The model 
focuses on the binary breakup of fragments formed by 
the primary cracking mechanism of the shell. Represent- 
ing the fragments by rectangles with a continuous mass 
distribution at a fixed aspect ratio r = a/b, the effect 
of the material dependent primary cracking mechanism 
on the shape of fragments can be taken into account by 
setting r w 1 for fragments of isotropic shape (egg), and 
r 3> 1 for the highly anisotropic needle-like pieces (glass). 
Based on Figs.n]and|21 these fragments are then assumed 
to undergo a sequential binary breakup process, where 
at each step of the hierarchy they break into two pieces 
of equal mass with a probability p < 1. Note that the 
fragments have 1 — p chance to keep their actual size. To 
capture the effect of out-of-plain deformations, we choose 
a side of a rectangle to break with a probability propor- 
tional to its length. Computer simulations of the model 
were performed starting from a continuous distribution of 
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FIG. 5: Log-log plot of (m) as a function of R g . Inset: dis- 
tribution / of the shape parameter S. 



fragment sizes varying the initial aspect ratio r to model 
different materials, while p was fixed. The hierarchical 
process was followed up to n = 30 generations resulting in 
~ 10 7 fragments in the final state, where the mass m, the 
radius of gyration R g and the shape parameter S of frag- 
ments were determined. It can be observed in Fig. [S] that 
similar to the experiments, for all values of r the average 
fragment mass exhibits a power law dependence on the 
radius of gyration (m) ~ Rg. Starting the process with 
isotropic shapes r ks 1, the fragments remain isotropic at 
all levels of the hierarchy implying a = 2 ± 0.05 as it was 
observed for the egg pieces. Modelling glass fragments by 
a high initial anisotropy r> 1, the crack mostly occurs 
along the same side of the rectangle lowering r, how- 
ever, when the two sides become comparable r fluctuates 
about one. Consequently, large fragments are character- 
ized by an exponent a significantly lower than 2, while 
for small pieces a crossover occurs to isotropic shape with 
a = 2. In the experiments we estimated the initial as- 
pect ratio of glass fragments from Fig. [21 to fall in the 
range of 15 < r < 35. Simulations with these aspect 



ratios r proved to provide values of the exponent a in 
the vicinity of 1.5. In Fig.[S]the results are presented for 
r = 20 where a = 1.5 ± 0.06 was obtained . When the 
fragments initially have an isotropic shape, the hierarchi- 
cal process gives rise to a rapidly decreasing distribution 
of the shape parameter f(S) over a narrow range as it 
was observed for egg pieces. Starting the simulation with 
elongated fragments (r = 20 in Fig. [5} the distribution 
/ shows a power law decay f(S) ~ S~P with an expo- 
nent (3 = 3.8 ± 0.3 very close to the experimental value 
of glass fragments. It is very important to notice that 
these results for the exponent a and for the distribution 
of fragment shapes are practically independent on p. 

Summarizing, based on high speed imaging tech- 
niques we have determined the fragmentation mechanism 
of closed shells: after the material-dependent primary 
cracking mechanism governed by the in-plane deforma- 
tion of the shell, a hierarchical secondary breakup pro- 
cess sets in due to out-of-plane deformations. Contrary 
to bulk systems, the shape of shell fragments shows large 
variations from completely isotropic to highly anisotropic 
fragments depending on the primary cracking mecha- 
nism. We pointed out that the anisotropic fragments 
have a sclf-affine character with a scaling exponent 6 — 
1/2. To give a quantitative characterization of fragment 
shapes we proposed a shape parameter the distribution 
of which was found to exhibit a power law decay. A 
hierarchical stochastic breakup model provided quanti- 
tative agreement with the experimental findings, which 
demonstrates the robustness of the scaling laws of frag- 
ment shapes in shell fragmentation. The results imply 
that the characterization of fragment shapes in breakup 
models of space debris production should be improved by 
a clear distinction of bulk and shell fragmentation and 
by using scaling laws with exponents depending on the 
cracking mechanism of the material. 
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